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Dear friends, 

there is so much to report, that I really don't 
know where to start. This time I have to do 
without a German version and I hope that my 
many German speaking friends will forgive me. 



The ACDCA-Summer Academy in Gos- 
ing, Lower Austria was a wonderful meeting of 
CAS users and -creators: David Stoutemyer and 
Theresa Shelby presented DERIVE Version 5, 
Frank Postel showed an interesting educational 
application of MuPAD and Danny Gremillion 
gave an impressive introduction to TI- 
Interactive. 65 delegates from 19 countries at- 
tended 35 lectures - a full program for three 
days, accomplished by a nostalgic train ride and 
a hike in the mountains. 

The picture shows Josef together with Carl 
Leinbach at the "End of the Trail". One or two 
cool drinks helped them to recover. 

Thank you all for coming and sharing your ex- 
perience. The papers will be presented on 
ACDCA's homepage very soon. We also plan to 
produce a CD with the submitted papers - de- 
pending on how many of the lecturers will de- 
liver their contributions. 

In this DNL you can find the "true Tidbits 
15". Unfortunately I repeated Johann Wiesen- 
bauer's Titbits #14 as #15 in the last DNL. 
Please study the new one carefully. It extends 
Richard Schorn's trip into Group Theory from 
DNL#33. (By the way, Richard took many pic- 
tures in Gosing and he promised to send some to 


put them also on the Summer Academy CD). As 
a reparation for my mistake I added two pages 
more of content and apologize once more for my 
mistake. 

Have you ever had the idea to calculate the 
"half derivative" of a function? David Halprin 
did. Follow his findings in the "Lighter Side of 
Operational Calculus". David also sent an amus- 
ing collection of "Epi-taphs/grams for Mathema- 
ticians of Desparate Fame". I hope to find some 
space in the next DNLs to present some pearls of 
his collection. 

"The News" of this summer were doubtless 
TLs take-over of DERIVE. Many messages full 
of fear, disappointment, but also full of hope and 
good ideas could be found in the news groups. I 
am happy to have the statements of two very 
competent people in the heart of this issue. Our 
good wishes for a prospering future accompanies 
the just married CAS-couple. 

In the next DNL you will find an interesting 
astronomic article dealing with the "3 bodies' 
problem", contributed by Ludwig Fritsch, 
Vienna. He uses DERIVE together with 
CYCLONE 98 to visualize this famous problem. 
Page 10 shows another CYCLONE 98 created 
object and announces its even more powerful 
successor. 

I'd like to draw your attention to the out- 
come of an European Union supported 
COMENIUS project. The book shelf contains an 
English and a German version [5]. It is very 
likely that you will know some of the project 
group's partners. 

You can see that it was a very busy summer for 
the DERIVIANS and Tl-ers. And I didn't even 
mention the Plymouth Conference - I couldn't 
attend this conference because we had to organ- 
ize the Summer Academy, I didn't mention the 
eclipse, I didn't mention 

Best regards until December, 

Josef 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & TI-92 User Group. It 
is published at least four times a year with 
a contents of 44 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE 
and the TI-92/89 as well as to create a 
group to discuss the possibilities of new 
methodical and didactical manners in 
teaching mathematics. 

As many of the DERIVE Users are also 
using the TI-92/89 the DNL tries to com- 
bine the applications of these modern tech- 
nologies. 


Preview: Contributions for the next issues 

3D-Geometry, Reichel, AUT 

Graphic Integration, Linear Programming, Various Projections a.o., Bohm, AUT 

A Utility file for complex dynamic systems, Lechner, AUT 

Examples for Statistics, Roeloffs, NL 

Quaternion Algebra, Sirota, RUS 

Various Training Programs for the Tl 

A critical comment on the "Delayed Assignation" :==, Kiimmel, GER 

Sand Dunes, River Meander and Elastica Halprin, AUS 

Type checking, Finite continued fractions, Winding number, Welke, GER 
Share Holder's Considerations using a CAS, Bohm, AUT 
Kaprekar's "Self numbers", GER 

Implicit Multivalue Bivariate Function 3D Plots, Biryukov, RUS 

LU-Decomposition, Morales, COL 

ODEs with Constant Coefficients, Fernandez, ARG 

Three Bodies Problem (DERIVE & Cyclone98), Fritsch, AUT 

Calendar Problems (Poblacion), ESP 

and 

Setif, FRA; Vermeylen, BEL; Leinbach, USA; Speck, NZL; Biryukow, RUS 
Wiesenbauer, AUT; Aue, GER; Koller, AUT, Ibrahim & Cordoba, ESP 
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e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to re- 
inforce the international touch of the DNL. 
It must be said, though, that non-English 
articles will be warmly welcomed none- 
theless. Your contributions will be edited 
but not assessed. By submitting articles 
the author gives his consent for reprinting 
it in the DNL. The more contributions you 
will send, the more lively and richer in 
contents the DERIVE & TI-92 Newsletter 
will be. 


Next issue: December 1999 

Deadline 15 November 1999 
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A Substitution Function for Variables 

by Stefan Welke, Spwelke@aol.com 
Abstract 

DERIVE has no built-in function to perform substitutions but offers the pos- 
sibility to substitute variables and subexpressions via menu commands. This article 
describes a function that performs substitution of variables. 


Introduction 

Consider the following problem: Given the equation of an algebraic curve y, e.g. y 2 - x = 0, find an 
equation for y in polar coordinates, i.e. r :=/( r). A handmade solution with DERIVE is: Enter the equa- 
tion, use the SUBSTITUTION - command to substitute x := r cos(t ) and y := r sin(t), then solve for r. 

Z 

#1: y - x 


2 

#2: Cr-SIN(t)) - r-COS(t) 

COS(t) 1 

r = 0, r = 

2 

SIN(t) J 

Comment 2015: Later DERIVE versions present the output of SOLVE in another form. Addi- 
tionally there is also the SOLUTIONS-function available: 


2 

SOLVEC(r-SINCt)) - r-CQS(t), r) = 

#3 : 


2 

SOLVE (( r - SI N(t)) - r-COS(t), r) 


2 

SOLUTIONSCfr ■ SIN ft)) - r-COSCt) 






COS(t) 

2 

SIN(t) 


i 

v r = 0 


J 


, r) 


COS ft) | 

0 , 

2 

SIN(t) I 


This way is particulary satisfactory for educational purposes because the student is aware of each 
step of calculation. But at a more advanced level one would like to automate this procedure. A solu- 
tion of this problem should be a function TO_POLAR (expr_) which returns an expression for /f r). 
The aim of this article is to provide a function SUBST which works like the SUBSTITUTION - pro- 
cedure for variables: 

2 2 2 
SUBSTCy - x, [x, y], [r-COS(t), r-SIN(t)]) = r -SIN(t) - r-COSOO 

The syntax is: SUBST (expr_, var_list, substitution_list) . In the expression denoted 
by expr_ each element of var_list is substituted by the corresponding element of substitu- 
tion list. The announced function TO POLAR now reads: 
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TO_PQLAR(exp r_ , var_> := RHS (SOLVE (SUBST(expr_, va r_, [r-COS(t), r-SIN(t>]> = 0, r» 


2 

T0_P0LAR(y - x, [x, y]) = 


0 , 


COS(t) 


2 


SIN(t) 

Comment: This is the TO_POLAR-function in DERIVE 6 syntax: 

to_polar(expr_, var_) :z SOLUTIONS (SUBST(expr_, var_, [r-COS(t), r-SIN(t)]), r) 


to_polar(y - x, [x, y]> = 


COS(t) 


r - 


v r — 0 


2 


SIN(t) 


Simple Substitution 

A list of expressions is evaluated from left to right. Thus the following simple construction works: 

r 2 1 r 2 l 

#4: Lx := y, y - x a x := xj = Ly a y - y, xj 

Note that we must clear variable x, because otherwise x keeps its value y for all future calculations. 
The second element of the RHS in #4 is what we really want and therefore we try the next command: 

f 2 2 1 Z 

#5: Lx := w s y - x, x := xj = y - x 

Z 

The reason for the failure is that only the second element of the LHS list of #5 is evaluated. A way 
around is to use a self-defined function which does nothing else but . This may look strange but it 
works because in user-defined functions the argument is evaluated first. 

PART(v_, n_) := v_ 

#6: m_ 

#7: 5IMPLE_SUBST(expr_ a var_, svar_) := PART ([ var_ := svar_, exp r_ a var_ := var_] , Z) 

3 3 

#3: 5IMPLE_SUBST(a- x + b, x s z ) = az + b 

This is exactly what we need as substitution of a single variable. 

Comment: #4 to #8 do not work in DERIVE 6. Use the implemented SUBST-function instead: 

3 3 

SUBSTCa-x + b, x a z ) = a-z + b 


Substitution of several variables 

We can use SIMPLE_SUBST repeatedly to substitute several variables in an expression. To do 
this we utilize the ITERATE-function. Consider the list v_ : = [ k , expr_] where k e N is a counter 
and expr_ is the expression on which we perform the substitution. We now iterate the following func- 
tion S over the argument v_. For brevity a_ denotes the variable list and b_ the substitution list. 


D-N-L#35 


Stefan Welke: A Substitution Function for Variables 


P 5 


#9: 


S(v_, a_ , b_) := Tl + v_ , SIMPLE_5UBST(v_ , a_ , b_ )1 
L 1 2 V_il V_il J 


Notice that the counter increases by one when S is applied to v_. We observe the effect of repeated 
application of S : 


#10: 

sd 


z 

X 


L [x, 

y, z], 1 

f 3 

Lw a 

a , 5 ] 

U 

#11: 

sd 

k 

6 

w 


L [x, 

y, z], 1 

r 3 

Lw a 

a , 5 ] 

U 

#1Z : 

sd 

k 

6 

w 

3 

- a z 

:], [ x , 

, y, z], 

[ 3 
Lw 

3 

, a , 5 

d> = 


6 3 

w - a z_ 


The function S2 uses ITERATE to produce the final result. The number of iterations equals the 
number of variables to substitute. S3 removes the counter: 

#13: SZ(expr_ a a_, b_) := ITERATE (S(w_ a a_ a b_) , w_, [1, expr_] , DIMENSION(a_)) 



z 

[331 

r e 31 

#14: 

S2Cx - yz, [x, y, z] , I 

Lw , a , 5j) = 

|_4, w - 5 a J 

#15: 

S3(expr_, var_, svar_) := 

PARTf52Cexpr_ 

, var_ a svar_) a 


^ I 

r 3 3 i 

6 3 

#16: 

S3Cx - y-z, [x, y, z], I 

Lw , a , 5_|) = 

w - 5- a 


Comment: #10 to #16 do not work in DERIVE 6 - and it is notnecessary. Use the imple- 
mented SUBST-function instead: 


SUBSTCx - y-z a [x, y, z] : 


l\ a 3 , 5 ]> = 


6 3 

w - 5- a 


Intermediate Variables 

Let us substitute x+ y for x and 2 + x for y in the product p := xy. The result should be 

P' := (x + y) (2 + x). 

But function S3 gives a different result: 

#17: S3(x-y a [x, y], [x 4 y 5 2 + x]) = 2-(x + 1) ■ (x + Z) 

The reason is that S3 works sequentially. First x is substituted, theny in both places: 
x y —> (x + y) y —> (x + (x + 2)) (x + 2) = 2 (x + 1) (x + 2) 

But only the second y should be replaced by x + 2. We can avoid this dilemma by the introduction 
of as many intermediate variables as that number of variables is that we have to substitute: 


#13: S3(x-y a [x, y, a s b] s [a, b, x 4 y, Z 4- x]) = (x 4 Z)-(x 4 y) 
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First the variables x and y are renamed. Doing so they cannot be confused with the variables x and y 
in the substituting expressions. Then the variables a and b are replaced by x + y resp. 2 + x. 

If we want to perform this procedure within a function we must ensure that the intermediate vari- 
ables are not accidentally used in that expression that is used for substitution. Since DERIVE does not 

have local variables a good way to prevent this difficulty is to use underscored variables , b , 

etc.: 

#19: v_ := [a b a c_, d , e_, f , g_, h ] 

The functions AUX_1 and AUX_2 enlarge the given variable and substitution lists by the correct 
number of intermediate variables from the list v_: 

AUX_l(u_ , v_) :z APPEND(u_ a VECTGR(v_ , i a 1, DIMENSION (u_))) 

#20: i 

#21: AUX_l([k , 1 , m , n] s v_) = [k , 1 a m , n a a , b , c , d ] 

AUX_2(u_ , v_) := APPEND (VECTOR (u_ s i a 1, DIMENSION(v_) ) , v_) 

#22 : i 

#23: AUX_2(v_ , [a, b, c]) = [a a b , c a a, b, c] 

These examples illustrate the use of AUX_1 and AUX_2. Now we are well prepared to define the func- 
tion SUB ST we were looking for: 

#24: SUB ST (expr_ , var_, svar_) := S3(expr_, AUX_l(var_ , v_) , AUX_2(v_ ? svar_)> 

#25: 5UBST(x-y, [x, y] a [x + y a 2 +- x]) = (x + 2) - (x -f y) 

2 2 2 
#26: SUBSTCy - x, [x, y] , [r-C0S(t) a r-SIN(t)]) = r -SIN(t) - r-COSCt) 

#25 and #26 show the expected results. The number of variables in v_ is a theoretical limitation to 
the application of SUBST. But v_ can be enlarged to math practically every purpose. Note that the 
number of elements of var_ must match the number of elements of svar_. 

Applications 

Our first application is the function TO_POLAR which was mentioned in the introduction. 


#27: T0_P0LAR (exp r_ , var_) := RHS(SOLVECSUBSTtexpr_ s var_ s [r-COSCt), r-SIN(t)]) = 0, r» 


2 

T0_PGLAR(y - x, [x a y]) 

#28: 


coset) I 

0 a 

2 

SIN(t) I 


2 2 

T0_P0LARCx - y +1, [x, y]) 

#29: 


1 


2 

l 7(1 - 2 ■ COS(t) ) 


1 


2 

7(1 - 2 ■ COS(t) ) J 


Comment: See the TO_POLAR-function in DERIVE 6 syntax on the next page: 
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f 


2 2 

to_polar(x - y +l,[x 1 y])= r 


v 


2 

5IGN(2 ■ CQS(t) - 1) 

v r 

2 

- 2- COSCt) ) 


2 

SIGN (2 ■ COSCt) - 1) 


2 

7C1 - 2- COSCt) ) , 


2 2 

to_polar(x +y -1, [x, y]) = (r =-lv r =1) 

f 1 i l 

to_polar(x- y -1, [ x , y ] ) = I r = - v r = 

l 7CsiNCt))-7CcosCt)) 7csiNCt))-7Ccosct)) , 


A second application is the inversion of an algebraic curve with respect to a circle. Formally inver- 
sion with respect to the unit circle consists of a substitution of variables in the defining polynomial of 
the algebraic curve: 


(x, y) 


y 


x + y 


+ y 


More mathematical background is presented in [1] und [2]. Consider the parabola y 2 - 2x = 0 as an 
example: 


#30: SUBST 


f 2 
y 


2x, [x, y] a 



3 2 2 

2x + 2-x-y -y 


2 2 2 
(x + y ) 


The numerator of #30 is the defining polynomial of the image curve under inversion. 


#31: 


NUMERATOR 


f 


V 


3 2 2 ^ 

2 x + 2 - x- y - y 


2 2 2 
Cx + y ) 


3 2 

= - 2 ■ x - y - (7 - x - 1) 


The following auxiliary function AUX computes the substitution list for the inversion with respect 
to a circle with radius r and center m . 


2 

r_ - Cvar_ - m_) 

#32: AUXCvar_, m_ 5 r_) := m_ + 

( var_ - m_) - ( var_ - m_) 


#33: AUXC[x 5 y] , [-2 a l] a 1) 


#34: 


x + 2 


2 2 

Lx i 4-x i y - 2 - y + 5 


- 2, 


2 2 
x +4-x + y - y + 4 


2 2 

x i 4-x i y - 2-y + Sj 


We combine all steps in one function. Note that we can use the function value as substitution list. 
Again the reason is that AUX is evaluated before SUBST. The FACTOR function is necessary to guar- 
antee that the result of SUBST has fractional form. 
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INVERSION(expr_, var_, m_, r_) := NUMERATOR ( FACTOR (SUBST(expr_, var_, AUX(var_, m_, r_))>) 
2 

INVERSION^ - 2-x, [x, y], [1, 0], 1) 

4 3 2 2 2 2 2 

- 2-x + 6-x - 2-x -(2-y + 3) + 2-x-(3-y f 1) - y -(2-y + 1) 



Remarks: My efforts to define a substitution function were 
motivated by the second example. When writing [1] I missed the 
possibility to use substitution of variables as part of a user- 
defined function. Mathematica and Maple offer such a function 
and I wondered if it was possible to write a function for DERIVE 
with the polynomial of the image curve under inversion as the 
function value. 

Further information on inversion of algebraic curves is acces- 
sible in [2]. 


Conclusions 

DERIVE r s strength lies in its restriction to a small set of really essential programming functions. This 
feature is particularly useful in education. But DERIVE r s programming functions are powerful enough 
to achieve a well working substitution function with only nine lines of programming code. So substitu- 
tion of variables can become a part of user-defined functions. Two additional lines of programming 
code provide the function INVERSION which manipulates polynomials in a nontrivial way. 


References 

[1] S. Welke, Inversion of Elementary Algebraic Curves with Respect to a Circle , Proceedings of the 
International DERIVE and TI-92 Conference, SchloB Birlinghoven, July 1996. ZKL Uni Munster 

[2] Xah Lee, http://www.best.com/~xah/SpecialPlaneCurves_dir/Inversion_dir/inversion.html 
(with many interesting links) 

Comment of the Editor: DERIVE 5 will provide a lot of additional programming tools. At the 
ACDCA Summer Academy David Stoutemyer and Theresa Shelby gave an exciting presen- 
tation of DERIVE 5 which will contain many new features. 

Comment for the revised version 2015: As mentioned above we have now the SUBST- 
command implemented. Stefan’s INVERSION-function works pretty well. I’d like to propose a 
small improvement. When plotting the inversion of the parabola given above you will receive 
an error message because you have to change to: 

4 3 2 2 2 2 2 

-2-x + 6-x - 2-x -(2-y +3)+2-x-(3-y + 1) - y -(2-y + 1) = 0 
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My proposal is: 

INVERSION (exp r_, var_, m_, r_) := NUMERATOR(FACTOR(SUBST (expr_, var_, 
AUX(var_, m_, r_)))) = 0 


Examples: 

Inversion of a straight line wrt to a circle [(1 ,0);1] gives a circle passing the origin and vice 
versa and 2 nd example is the inversion of an ellipse wrt to the circle [(1 ,0);2] and vice versa. 

INVERSIONS- x + y - 1, [x, y], [1, 0], 1) 

Z Z 

x + y + y = 1 

Z Z 

INVERSIONS + y + y - 1, [x, y] , [1, 0], 1) 

Z- x + y = 1 



Z Z 

(x - 1) +■ y =4 

Z Z 

INVERSIONCZ-x + 10 -y - 15, [x, y] , [1, 0], Z) 

4 3 Z Z Z 4 Z 

13 x - 68- x + Z-x ■ (13 ■ y + 47) - 4 x (17 y + 9) + 13-y - 118-y = 3 
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Morris WaehltEs M Ten Section Method” 

In July I had a telephone talk with Morris who wanted to create a DERIVE - program to demonstrate 
an algorithm for finding a root of an equation numerically similar to the Bisection-Method, but using a 
"Ten Section". Some days after I had a FAX, which described accurately what he had in mind: 


As per our telephone conversation, here, on paper, is what I was thinking about. To create a DERIVE 
’’program” to find the square root of 5 (or any other number) by the TEN SECTION METHOD. 


NUMBER 

NUMBER belongs 
to VECTOR 

SINCE (test) 

1 

no 

l 2 < 5 and 2 2 < 5 

2 

yes 

2 2 < 5 and 3 2 > 5 

2.0 

no 

2.0 2 < 5 and 2.1 2 <5 

2.1 

no 

2.1 2 < 5 and 2.2 2 < 5 

2.2 

yes 

2.2 2 < 5 and 2.3 2 > 5 

2.20 

no 

2.20 2 < 5 and 2.2 1 2 < 5 

2.21 

no 

2.21 2 < 5 and 2.22 2 < 5 

2.22 

no 

2.22 2 < 5 and 2.23 2 < 5 

2.23 

yes 

2.23 2 < 5 and 2.24 2 > 5 

2.230 

no 


2.231 

no 


.... and so on .... 




The vector in the PROCESS of GROWING would look like this: 

[ 2 , 2 . 2 , 2 . 23 , 2 . 236 , 2 . 2360 , ] 

It would STOP GROWING at some desired degree of precision. 


If one wanted to find V59 , perhaps one could give the ’’program” 7 as an initial value to speed the 
process up a bit. 

Finding V5 is like solving the equation x 2 - 5 = 0 for its positive root. With some modifications to the 
’’program” could one use it to find a root of some less simple equation f(x) = 0 by the TEN SECTION 
METHOD, perhaps by giving it an initial value?? A number, say 3, would belong to the solution vec- 
tor if f(3) < 0 but f(4) > 0 i.e. 


NUMBER 

NUMBER belongs 
to VECTOR 

SINCE (test) 

3 

yes 

f(3) < 0 and f(4) > 0 

3.0 

no 

f(3) < 0 and f(3.1) < 0 

3.1 

yes 

f(3.1) < 0 and f(3.2) > 0 

3.10 

no 

f(3.10) < 0 and f(3.11) < 0 

3.11 

yes 

f(3.1 1) < 0 and f(3.12) > 0 

3.110 

yes 

f(3.1 10) < 0 and 3.1 1 1 >0 
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The solution vector in the PROCESS of GROWING would look like this: 

[3, 3.1, 3.11, 3, 110, ] 

It would STOP GROWING at some desired degree of precision. 

I appreciate that the TEN SECTION METHOD is NOT a very efficient way to do anything. BUT I am 
regarding it more as an EXERCISE in programming with DERIVE. As a TEST, if you like, of what 
one CAN and CAN NOT do with DERIVE functional programming. 

For instance, is it possible to have a DERIVE Program PAUSE and accept user input??? 

(Hi Al, Dave and Theresa, it is to you to answer this question!) 

If you think that all this will turn out to be just a time consuming bother, PLEASE just forget about it. 
Best regards, Morris Waehlti, Geneva, Switzerland. 


I didn't forget and tried my best, but I was not as successful as with the "Look and Say Se- 
quence" (page 34): 

-k_ 

TENSEC(u, x, n) :z VECTOR(MAX(SELECT(u < 0, x, VECTOR(0 + n_-10 , n_, 0, 

k_ + 1 

10 ))), k_, 0, n) 

2 

TENSEC(x - 5, x, 4) 

[2, 2.2, 2.23, 2.236, 2.236] 
in 1999 80 sec, now only 6 sec 
2 

TENSEC(x - 5, x, 5) 

[2, 2.2, 2.23, 2.236, 2.236, 2.23606] 

in 1999: I couldn't wait longer!, now: needs 133 sec 


As Johann Wiesenbauer - my man for all problems connected with programming in 
DERIVE - was in Hawaii to work together with the Soft Warehousers to improve 
DERIVE 5 I wrote to one of his competitors in programming skills Terence Etchells 
and posted Morris' TEN SECTION PROBLEM. Very soon I received the following 
file: 
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"Ten Section Method" 

"Written by Terence Etchells, Liverpool John Moores University, UK" 

"For Morris Waehlti and the DUG" 

Precis ionDigits : =10 

Notation : ^Decimal 

SIGNS (w) : =IF (w=0, 1, SIGN (w) ) 

TEST_FOR_ROOT (f,x,a,b) :=IF (SIGNS (ITERATE (f,x,a, 1) ) =SIGNS (ITERATE (f,x,b, 1) ) , 
0,1,0) 

FIND_NEXT_SECTION ( f , x, s_, i) : ^ITERATE (IF ( TEST_FOR_ROOT (f ,x, a, a+i) =1, a, a+i, 
a+i) , a, s ) 

TEN_SECTION ( f , x, s_int, s_val,n) : ^ITERATES ( [i/10, 

FIND_NEXT_SECTION (f ,x, s_, i) ] , [i, s_] , [s_int, s_val] ,n) ' SUB 2 

f is the equation f(x) = 0; x is the variable; s int is the initial interval size, s val is the starting point; 
n is the number of times you wish to perform the decimal search, 

e.g. x 2 - 5=0, starting at x = 1 with a starting interval of size 1,10 times. 

TEN_SECTION (x A 2-5 , x, 1 , 1 , 10 ) 

; Simp (#13) 

[1,2,2.2,2.23,2.236,2.236,2.23606,2.236067,2.2360679,2.23606797, 

2.236067977] 

N.b. you must use a high precision level if n is large! ! ! 

E.g. x 2 - 5=0, starting at x = -20 with a starting interval of size 10, 10 times. 

TEN_SECTION (x A 2 -5 , x, 1 0 , -2 0 , 1 0 ) 

; Simp (#17) 

[-20,-10,-3,-2.3,-2.24,-2.237,-2.2361,-2.23607,-2.236068,-2.236068, 

-2.23606798] 

The interval size is divided by 10 for each iteration; this can be changed by editing the 
TEN_SECTION ( ) function where the i/ 1 0 appears. 


This nice object consisting of 27 spheres - 
it is really only ONE surface - was pre- 
sented by David Sjostrand at the 
ACDCA-Summer Academy 1999. 

He used DERIVE to edit and simplify the 
expression below and imported this expres- 
sion to Cyclone98 to plot this box full of 
spheres. 

Ill 2 2 2 

n n n ((x - 2*a) + (y - 2*b) + (s - 2*c) - 1) = 0 

c=-l b=-l a=-l 

There is a wonderful new version of Cyclone99 - DPGraph2000 (see page 33): 

http : //www . yds a . s e /DP Graph2 000/DP Graph2 000 /DP Graph2 000 . html 
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Peter Witthinrich, Liibeck, Germany Pwitthinrich@t-online . de 

Dear Josef, 

Alfonso Poblacion's algorithm shows another way to compute pi, but what is the basics, the theoretical 
background? 

I didn't intend to show a high speed algorithm but the spigot algorithm found by Rabinowitz/Wagon. 
And Poblacion's algorithm is not a spigot algorithm, which produces the decimal places digits by dig- 
its (so the digits ’drop’ one after the other). 

Another important difference lies in the fact that the spigot algorithm uses only integer arithmetic, i.e. 
there are no rounding errors apart from integer overflow (which can be avoided by suitable routines or 
exception handling). So the spigot algorithm uses low level mathematics, but I think the PIAPPROX 
algorithm uses Derive' s internal floating point routines. They are good I know, they were already 
good in old muMath years ago. Anyway the TI-92 is much slower and lacks memory for computing 
large numbers. 

My proposal: Alfonso Poblacion describes us the theoretical background of his mentioned algorithm. 
There are several algorithms computing pi: I would like to learn about another one. 

Best wishes 

Viele GriiBe aus Liibeck 


Alfonso Poblacion, Valladolid, Spain alfonso@ gauss .mat . eup . uva . es 

Dear Josef and Peter, 

By the days I read DNL #33,1 was working with other colleagues on the writing of a Derive 
Secondary level book. In this book we include several miscellaneous chapters combining maths, his- 
tory and courious things. One of them is about squaring the circle and the calculation of PI. So I read 
some articles and as I found interesting (because of the speed) algorithms to compute PI, I sent you 
one of them for the DNL relating this with the appearance of the Spigot algorithm. But my aim was 
not to say something like "This is a better algoritn because is faster". I know both are different, so 
please Peter, forgive me if I caused you that bad impression, surely deriving from my poor English. 

The algorithn (that Josef gently named PIPOB, but of course is not mine) was taken from an 
article in the Spanish version of Scientific American titled "Ramanujan and PI" that appeared in Spain 
in the issue of April 1988, whose authors were Jonathan M. Borwein and Peter B. Borwein. They ex- 
plained on that article that they deduced the algorithm modifying another one from Gauss that was 
later modified by Richard Brent and Eugene Salamin (called Gauss-Brent-Salamin algorithm). After 
the Borwein & Borwein's changes, they said that was used by Yasumasa Kanada in 1983 to compute 
thousands of digits of PI (As you probably know, Kanada is an expert and recordman about calculat- 
ing Pi digits). So the explanations I can give you are the ones from that article that also includes an- 
other two iterative algorithms. 

They explained that these algorithms are based on modular equations. A modular function 
L(q) is one that can relate by an algebraical expression called modular equation involving the same 
function expressed by the same variable q raised up to an integer power L(q A p). The integer p deter- 
mines the "order" of the modular equation. For example (I will write mathematical expressions like if 
they would be written with DERIVE, but do not try to translate them to DERIVE; I only use it like a 
word processor) 

L(q)=16q PRODUCT ((( l + q A ( 2n ))/( l + q A ( 2n-l ))) A 8 , n, 1, inf) 

Its associated equation is of 7 th order and it relates L(q) with L(q A 7): 

( L ( q) L ( q A 7 ) ) A ( 1/8 ) + ( ( 1-L (q) ) (1-L(q"7) ) ) A (1/8) =1 
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The solutions of the modular equation that satisfy added conditions are called singular values. One of 
these values is deriving from a sequence of values K sub p (I don not know how to use this editor to 
write subindices), where 

K sub p = sqrt (L (e A (-pi sqrt (p) ) ) ) 


These values has the property that -2/sqrt(p) log(K sub p 14 ) coincide with a lot of digits of PI. The 
number of common digits increases with p. 

Ramanujan used this general procedure to build some series related with PI. On the algorithm 
described in the DNL, the first step (the calculation of y(n)) belongs to the calculation of one of the 
terms of a sequence of singular values that is obtaining when you solve a modular equation of a suit- 
able order (the authors did not say which). The 2 nd step (calculation of a(n)) belongs to the logarithm 
of the singular value. 


This is more or less what the authors say about the background of the iterative algorithms in 
the referred article. I hope this would be what you asked me for. 


Have excellent holidays. 
Alfonso J. Poblacion. 


DtL: A. DERIVE to LaTeX-Translator 

Gary Bishop, bishop@cs.unc.edu 
August 1999 


1 What is it? 

DtL is a little program I wrote to address the problem of translating equations from Derive into the 
form required by LaTeX. Manually translating anything more than trivial expressions is error prone 
and frustrating. With DtL I can work out an equation in Derive and immediately paste it into my La- 
TeX document to get publication quality formatting. 

2 How do I use it? 

In Derive for Windows, select an expression or sub-expression and copy it to the clipboard with "Copy 
Expressions' from the "Edit' menu or by holding the Control key and typing C. Next run DtL.exe; a 
shortcut on the desktop, or a "shortcut key' is great for this. DtL will replace the Derive expression on 
the clipboard with the equivalent input for LaTeX. Finally, switch to the editor you use for entering 
LaTeX and paste the result. 

You can also use DtL with the " .mth' files produced by Derive. Run DtL with the command line: 
DtL.exe -f -F < myfile.mth > myfile.tex to convert all the expressions in myfile.mth into LaTeX for- 
mat. 

3 How does it work? 

When you select an equation in Derive for Windows, the text of the equation is placed on the clip- 
board in two different formats. The format you have access to in most programs is formatted text us- 
ing the DfW printer font. You can paste this form into a document to get a representation just like you 
saw in the Derive window but the quality of the result is well below what you get with a serious type- 
setter like LaTeX. 

The other format, "DfW Expression' is the one operated on by DtL. When you run DtL it reads the 
"DfW Expression' from the clipboard and translates it into an equivalent LaTeX representation. It then 
exits leaving the input for LaTeX on the clipboard. Then you can paste this text into whatever text 
editor you use. I prefer Emacs with the Auctex package but any editor that supports pasting text from 
other applications should work. 
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I used flex and bison from the gnu family of compiler building tools to do the lexical scanning and 
parsing of the input. These are the free software equivalents of lex and yacc from Unix. The remainder 
of the code is pretty straight-forward C++. It would be easy to modify DtL to produce input for TeX or 
Matlab. 

4 How do I get DtL? 

Although I wrote DtL to meet my own needs, I have made the source and executable available at 
http://www.cs.unc.edu/~gb/DtL.html. You can download either a zip file or a tar.gz file. Download 
the format that you prefer, unpack it with a program like WinZip or pkunzip and put the executable 
DtL.exe somewhere on your path. I put a shortcut on my desktop so that I can easily invoke DtL 
whenever I need it. 

5 Where can I get the other software mentioned here? 

Emacs, the world's greatest and most programmable editor, can be downloaded from 
http://www.cs.washington.edu/homes/voelker/ntemacs.html. The auctex package customizes emacs 
for editing TeX and LaTeX; you can download it from http://sunsite.auc.dk/auctex/. The MikTeX 
package provides the LaTeX processor and other programs for previewing and printing; you can 
download it from http://www.snafu.de/~cschenk/miktex/index.html. All of the gnu tools I used were 
part of the cygwin package of Unix tools for Win32 available from http://sourceware.cygnus.com. 

6 Who supports DtL? 

Not me. If you find bugs in DtL, please send me email. I am likely to fix it but I won't promise to fix 
any problem. Feel free to fix it yourself and let me know what you changed; that is why I included the 
source. As I said earlier, extending or modifying DtL to change the output style or language should be 
easy for anyone who knows how to use the tools. If you have never used flex (lex) or bison (yacc), I 
highly recommend you learn to use them. These powerful tools make translators such as DtL a snap to 
program. 


File translated from T E X by TtH, version 2.30. 
On 03 Aug 1999, 15:36. 


Gerhard Bardach 


g . bardach0gmx . at 


Does someone know a function in Derive to create a chain of generalized eigenvectors? 
I need it for the creation of invariant subspaces for a matrix. 


Sebastiano Cappuccio, Italy p . cappuccio@fo . nettuno . it 

Dear Joseph, 
please, see the URL: 

http://people.ciram.unibo.it/~barozzi/Net Elenco schede.html 


Here everybody can download a whole book of exercises of Calculus 1 using Derive in PDF format. 
Unfortunately files are in Italian only. 

I look forward to meeting you in Goesing next August. 

Best regards 
Sebastiano 


The URL given above is valid. See also http://people.ciram.unibo.it/-barozzi/MI2/PDF/ 
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You are friendly invited to visit the various home pages of T -Europe 
institutions. They provide a lot of CAS-materials and informations and 
you can also find many interesting links to related sites, (all from 
2015) 

United Kingdom and Ireland 

http : / / www . tcubed . org . uk/home/ 

Finland: 

https : //education . ti . com/f i/ suomi /nonproducts ingle/opetta j a_koulutus 

Sweden: 

https : //education . ti . com/ sv/ sverige /larare/larar f ortbildning/kostnadsf ria- 
larartraf f ar 

Germany: 

http : // www . t3deutschland . de/ 

Switzerland (German & French) 

http : //www . t 3 schweiz . ch/ 
http : //www . t3 suisse. ch/ 

Austria 

http : // www . t3oesterreich . at/ 

France: 

https : / / education . ti . com/ f r/ trance/ formations/ t3f ranee home 

Portugal: 

https : //education . ti . com/pt/portugal /professor /professor 
http : //www . apm . pt/apm/t3/t3 . htm 

Italy: 

https : // education . ti . com/ it/italia/professori/ aggiornamento 


T3 Worldwide 

education.ti.com/en/us/pd/community/t3-ww 



T 3 t U B D P E 


United States and Canada: 

https : //education . ti . com/en/us/pd 

Belgium (Flemish and Walloon): 

http : / / www . t3vlaanderen . be/ 
http : // www . t3wallonie . be/ 

Holland 

http : //www . t3nederland . nl/home/ 

Denmark 

https : //education . ti . com/da/danmark/nonproductsingle/laerere t3 

Norway 

https : // education . ti . com/ no / no rge/ nonproduct single/ it 3 

Australia and New Zealand 

https : // education . ti . com/ en-GB/ aus-nz/pd/ online -learning 
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Where B2n is equal to the Derive XM definition of BERNOULLI(2n) in the file NUMBER.MTH 
and is the nth term in some tables of Bernoulli numbers. 

d' a .sinx = sunx , D^.cosx = cusx , D^-tanx = tunx 

D^.cosecx = cosucx , D^secx = sucx , • cotan x = cotun x 

1 _ 1 _ sun x _ 

COSUC2 x 5 SUC2 X > tun2 x 

sun x cus x cus x 

1 _ 1 _ sunh x _ 

; cosuch 2 x j : such 2 x ? : tunh 2 x 

sunh x cush x cush x 

D^.ln (1 + x) = lun (1 + x) 



TANx TAN x and TUN x 


This shows that the second semi-derivative of sin x converts it to cos x, while the second semi- 
derivative of cos x looks different from what we might expect, so we substitute (r+1) for r and we have the 
correct expression, save for the exponent of '-T effecting a change in its sign, as expected. 
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TANH x = SINH x / COSH x 


SUNH x / CUSH x 




SINH x, SUNH x COSH x 
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SUNH x and CUSH x 


LN(1 + x), LUN(1 + x) and l/(l+x) 


Now to a logarithmic series: 


ln(l + x) = £ 


^-(-l) r -x r 


D^.ln(l + x) = y ' ( ~ 1)r - r! - X ^ 

U r.(r-^)! 


D*.D*.fa(l + x) = ° I' (-D' x- 1 = Q.l 

f7 r.(r-l)! j-f 1 ^ v 


Now to some tests for distributivity, associativity 
and commutativity: 


n K r_ r! 'X r '^_ 

D 2 x 7— --F: 


(r-/ 2 )! 

(r-X)-r!-x r " LS _ r!. x r " L5 
(r-X)! (r-1.5)! 


D.x r = r.x rl 


1 / r tr-l)!.x r ' 1 ' 5 r' x r ' 15 

.(D , x r ) = 1 — = — 

(r-1.5)! (r-1.5)! 

n i .5 r__ r!-X r_1 ' 5 

D x_ (Ud)i 

Now to try to find a chain rule for a product of two functiohs:D^ • D = D . W 2 = D 1 5 


m n _ m+n 
X X — X 


e.g. D .sinx = D /2 *cosx etc. 


T\ / m n\ m -1 n i m n -1 — m -1 n -1 _i_ 

v\x -x j-mx x + x nx -x x (mx + nx) 


= x m+n " 2 (m + n).x = (m + n). x m+n l Q.E.D. 


D 2 ■ x 


(m + n)! x m+p ' > 

(m + n-K)! 


,/ ,/ x n .m!. x m ^ x m .n!. x n '^! 

x D /2 x + x D /2 x ; — 7 — + — ; 7777 — 


(m-X)! 


(n-X)! 


_ m+n-y 
- X 7 


(m -X)! (n -X)! 


^ D /2 ■ x 
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^(1/2)! (3/2)! (5/2)! (7/2)! 

E(k) = Vi.(sun x -i.cus x) 


E(- ix ) = Vx vi. 


• 2 3 

X lx X 


v ; U 1/2 ) ! ( 3/2 ) ! ( 5/2 ) ! ( 7/2 ) ! ) 

E(-ix) = Vi.(i .sun x - cus x) 
i.E(ix) + E ( -ix) = 2iVi.sunx 
i • E ( ix) - E ( - ix) = 2 Vi .cus x 

cus 2 x + sun 2 x = - E(ix)-E(-ix) , cus 2 x-sun 2 x = )^i (e, 2 x) - E 2 ix) ) 



SIN x and SIN 2x 


d( e mx ) 


SUN x and 2 SUN x CUS x 


= m. e 


n 2 mx _ d 2 (e mx ) _ 2 mx 

D e : — ; m e 

dx 

Hence 


W / 2 mx — % y ( mx r 2 
D ' e " m if (r-X)l 


D K • D 1 ' 2 . e mx = m . Y (mx } = m . e mx Q.E.D. 

£ (r-1)! 

2 sun x .cus x = \ (E 2 ix) + E 2 ix) ) 


Now to demonstrate graphically one-third 
derivatives. 



SIN x, D A (l/3) SIN x, D A (2/3) SIN x and COS x 
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See a short section of David's HEAVISID.MTH 


#1: File HEAVISID. MTH , (c) 1998 by David Halprin, PO Box 460, Nth. Balwyn 

#2: Melbourne, Victoria, Australia. A file to investigate fractional derivatives. 

#3 : Notation := Rational 


n 2-n + 1 
10 C-l) -x 

#4: I 

n=0 (2-n + 1)1 

#5: Simplify to receive the series for sin x 

n 2-n 

10 C-l) ■ x 

#6: I 

n=0 (2-n)! 

#7: Simplify to receive the series for cos x 

n! n - 1/2 



#9: This is the operator to obtain the semi-derivative of x A n. 

#10: Now to apply it to the general terms in the trig series. 

n 2-n + 1/2 
10 C-l) -x 

I 



#12: 

#13: 

#14: 


#15: 


Simplify to obtain the semi-derivative of sin x = SUN x 
2-n + 1 

10 x 

I 

n=0 C2-n + 1)1 

This is the Taylor series for sinh x to 11 terms. 


2-n + 1/2 

10 x 



#16: Simplify to obtain the semi-derivative of sinh x = SUNH x 


BIBLIOGRAPHY 

Heaviside, Oli. Proceedings, Royal Society of London, 

Pt.I Feb 2 1893 Pp 504/529, Pt.II June 15 1893 Pp 105/143 

Poritsky, H. American Mathematical Monthly June 1936 Pp 331/344 
Peacock, G. British Association Report 1833 Pp 211/221 & 241/249 

David Halprin, P.O.Box 460 
North Balwyn, Victoria 3104, Australia 
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It was a nice day in August. I opened my electronic mailbox and found: 

On Tue, 10 Aug 1999, Terence Etchells wrote: 

Hello Fellow Derivers 

I strongly suggest that you go to the site 
http : //www . ti . com/ calc /docs/ news/ cg-504 . htm 
for some rather Earth shattering news. 

Texas Instruments Acquires DERIVE 

Many angry, sad, friendly, disappointed, satisfied and expectant messages followed. Let's read the 
opinion of two important and competent people from the two companies: 

To all DERIVE Users: 

On 6 August 1999 Soft Warehouse, Inc. and all its assets, including DERIVE, were acquired by Texas 
Instruments Incorporated. There has been much concern expressed on this NewsGroup and elsewhere 
about the impact this will have on the future of DERIVE. As one of the principal authors of the sys- 
tem, I would like to respond to these concerns. The views expressed are my own, not TPs. However, 
note that I am working with TI as an independent consultant. 

At the risk of over anthropomorphizing DERIVE, I make the following metaphor: More than one per- 
son has told me that DERIVE has many of the traits of a woman (intuitive, persistent, eager to help, 
etc.). I have been working on DERIVE and her predecessor muMATH for 20+ years now. And like 
any 20 year old, she is ready to go out into the world on her own. Like any father would be, I am 
proud and elated that she has come of age, but also deeply concerned about her future. 

The following facts have convinced me that passing the DERIVE torch to TI is the right thing to do: 
Soft Warehouse, mostly through the efforts of my partner David Stoutemyer, has had a close working 
relationship with the Educational and Productivity Solutions (E&PS) division (formerly the Calculator 
division) of TI for more than seven years. 

During that time E&PS has successfully developed and marketed products (the TI-92 and TI-89 calcu- 
lators) based in large part on computer algebra (CA) software developed in cooperation with Soft 
Warehouse. 

The investment TI made to acquire Soft Warehouse provides clear evidence that TI appreciates the 
value and potential of CA in general and DERIVE in particular. 

TI faces strong competition in the graphing calculator market. To remain competitive in this market, 
TI will need to provide math and science educators with an integrated set of math tools (both hardware 
and software) available on a broad range of platforms. 

Writing and debugging completely new software for each product would be prohibitively expensive, 
difficult to support, and not provide the required integration between products. Instead, it makes much 
more sense to write an interface specially designed for each product, but to base all the products on a 
single math engine. TI appears to want to base its products on the math engine that powers DERIVE. 

Theresa Shelby designed (in collaboration with me) and implemented the DERIVE for Windows inter- 
face (both version 4 and the forthcoming version 5). Theresa and David Stoutemyer are now employ- 
ees of TI. I am a consultant working with TI. Presumably we were hired by TI for our expertise in 
developing compact, efficient CA systems, as well as for our experience in marketing mathematical 
software over the last 20+ years. Therefore, we have a unique opportunity to significantly influence 
the direction E&PS takes in the future. 

Our agreement with TI makes it clear that the top priority for Theresa and me is to complete, thor- 
oughly test, and release version 5 of DERIVE for Windows as soon as possible. Thereafter I will con- 
tinue to work on the math engine that powers DERIVE. 
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Freed from the responsibilities of running a company, I can devote more of my professional time and 
creative efforts to extending the mathematical and programming powers of DERIVE. 


TI has the resources, personnel, and know-how required to effectively promote DERIVE in the current 
highly competitive global software market in a way that Soft Warehouse could never hope to do. 


In addition to its own considerable programming talent, TI has the stature to attract first-rate pro- 
grammers and mathematicians to enhance DERIVE's performance and reliability, and to extend the 
range of its capabilities. Also authors will be more likely to write and publishers more likely to pub- 
lish math and science books based on DERIVE if TFs name is on the product. 


Several users have expressed concern that DERIVE will be transformed from an efficient and compact 
CA system into a behemoth much like our competitors'. E&PS pioneered development of the graphing 
calculator, and more recently (in cooperation with Soft Warehouse) has produced calculators having 
symbolic math capabilities. The relatively modest computing power and memory size of such port- 
able platforms gives E&PS a vested interest in preserving DERIVE's compact size and high perform- 
ance. Also, the calculator paradigm that DERIVE uses fits in well with that used on products designed 
for math and science education. 


muLISP (the LISP system in which DERIVE is written) and DERIVE (and its predecessor muMATH) 
have steadily evolved over the past 20+ years. The mistakes and successes I made in the process have 
taught me a great deal about how to implement a computer algebra system. I am not going to live 
forever. I am determined to pass this knowledge on to others more capable than me so that DERIVE 
and/or her successors will continue to be used and improved. Thus, as soon as version 5 is released, 
one of my primary consulting duties will be to transfer this knowledge to software engineers at E&PS. 
I and the other authors greatly appreciate your loyalty and enthusiasm for DERIVE, and your concern 
about her future. But, now it's time to get back to work on version 5! 

Aloha, 

Albert D. Rich 
Co-author of DERIVE 


P.S. We recently received the following email from my friend and colleague, Julio Valella of Texas Instru- 
ments: 


Al, David, and Theresa: 

I am delighted to be assigned as the Texas Instruments DERIVE Product Manager. 

In my new role I have a monumental challenge, to continue your long established excellence in devel- 
opment and support of the Derive family of products. And to help you complete Derive 5.0 and ship, 
promote, and support it worldwide with our renown Tl-Cares™ Educator Support services. 

We are a mere few days into this challenge! We are learning, organizing, training our staff, and mar- 
shalling resources required to maintain uninterrupted delivery and support of DfW 4.1 1 while we pre- 
pare to launch and support DfW5.0 worldwide. Some of the forthcoming Tl-Cares ™ services Derive 
enthusiasts should enjoy in the near future include: 

+ Workshop Loan Program - Texas Instruments free loan service to educators and professionals who 
wish to conduct training, academic presentations, teacher workshops, and on site institutional dem- 
onstrations. 

+ Conference Exhibits - Texas Instruments exhibits at most math and science conferences worldwide. 
Beginning this November Derive will be displayed with opportunity for hands on demonstrations, 
literature, applications and support services information. 

+ World Wide Web information dedicated about Derive will continue at http://www.derive.com/ 
Soon global exposure and support for Derive expands at http://www.ti.com/calc/docs/software.htm 
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As the Team Leader of TI-92 development my emotional connections with Derive and its develop- 
ment team run deep. Now I am eager to get to know and work with Derive customers and enthusiasts 
worldwide. 


With a little help from all our friends I'm confident we will meet our challenge and measure up to the 
high expectations you have conditioned in Derive users! 

Cheers, 

Julio Valella 

Derive Product Manager 

Texas Instruments Educational & Productivity Solutions 

jvalella@ti.com 

http://www.ti.com/calc/ 


This wrote Heinz Rainer Geyerfrom Germany: 

Hey Albert, 

A wonderful letter. 

it's all about math, we know, but your comparison really got something! 

So, what is our (DERIVERS) play: I thing we all fear to loose that beautiful woman student, we (they 
all) all in love with, to that tall, rich and smart guy, who soon will leave her alone, broken and ruined. 

But sometimes that guy appears to be her real big love, and when they got married the future will start 
in their baby with even greater capacities: the cleverness of the grandfather combined with the 

smartness of the father. 

Back to reality: a few days ago at ACDCA conference in Goesing Danny Gremillion of TI presented 
some features of that new TI InterActive, which has no CAS in Version 1 and some other important 
features are also missing, but it looks like that tough guy. 

And DERIVE 5 appears dressed up, so he has to look at her. 

In our days education always has to have a look on fashion. Pupils got to have new designs like work- 
sheets with free shifting entries and one-click-connection to embedded objects, it also should have 
internet facilities to present your work on schools homepage and all that snick-snack. 

But one day the both (TI-InterActive and DERIVE) will have their baby, which usually get the name 
of that grandfather. 


Name it DERIVE InterActive 

So, wipe your tears away and take a drink with the grandfather in law. 

Cheers 

Rainer 


The last words again to the person, who started and encouraged the discussion on the inter- 
net, Terence Etchells: 

Well said Sir ! ! ! 

Cheers 


Terence 
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ON THE SOLUTION OF A LINEAR DIFFERENTIAL 
EQUATION OF THE ORDER n WITH CONSTANT 

COEFFICIENTS (3) 

Emilio Defez Candel , Vicente Soler Basauri , Departamento de Matematica Aplicada 
Universidad Politecnica de Valencia, P. O. BOX 22.012. Valencia. Spain 

4 Complete linear differential equation of the nth order with 
constant coefficients. 

Function DIFERENCIAL_HOMOGENEA_COEFICIENTES_CONSTANTES (v, x) 

can also be used to calculate solutions of the complete linear differential equation of the nth order with 
constant 

coefficients. 

a n y n) (x)+a n _ x y n ~ l \x)+ . . . + a x y\x)+ a 0 y(x) = f{x) . (20) 

Let us suppose that y ^(x) is a particular solution of (20). Then the change y — Z + y p transforms 

(20) into a type (6) homogeneous differential equation. Let us suppose that is the general 

solution of the homogeneous equation associated to (20) given by 

a n y n) (x)+a n _ x y w_1) (x)+ ... + a x y\x)+ tf 0 y(x) = 0 (21) 

It is then inferred that the general solution of equation (20) is given by 

>'(*) = y H U) + y P (*) ( 22 ) 

This solution contains the n arbitrary constants that must appear in the general solution of an differen- 
tial equation of the nth order, since they will be contained in the expression of y H {x). Based on the 
theorem of existence and uniqueness of solutions 2.1 it can be demonstrated that (22) is actually the 
general solution of (20). Thus, the problem of solving (20) is reduced to the problem of finding a 
y p {x) particular integral of (20). We are going to see two methods, the method of undetermined coef- 
ficients and the method of variation of constants, which will enable us to obtain this y p {x) particular 
solution of (20). 

4.1 Method of undetermined coefficients. 

This method consists of writing a similar expression to that of f(x) as an y p {x) possible particular 

solution, but with undetermined coefficients. These coefficients will be determined in such a way that 
when y is substituted by y (x)m (20), both members of the equality are identical. The following 

rules enable the each form. y p {x ) ac cording to the function f(x). 

1 . If fix) = P m (x) , where P m (x) is an m-degree polynomial in x, we will state 

+A x x + ... + A m x m ^x k (23) 


k being the multiplicity order of 0 as a root of the characteristic polynomial associated to equa- 
tion (21) and the coefficients {A 0 ,A l ,...,A m j are the constants to be determined. 
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2. If f(x) = e ax P m (x) , where Pjx) is an m-degree polynomial in x and aG ^, we will state 


y P { x )- eax (a + A x+ "• + A m xJn ^) xk 


(24) 


k being the multiplicity order of a as a root of the characteristic polynomial associated to 
equation (21) and the coefficients {A 0 , A x , . . . , A m } are the constants to be determined. 


P v (x)sin(p x) 

P c (x) cos(p x) , where in each one of the cases P s (x) and P c (x) 

P s ( x )sin(p x) + P c (x) cos(p x) 

c according to 
max {?,<:} 


3. If /(*)= 


are polynomials in x in s and c degrees and (3 e R , let m be m = 


which the one corresponding to /(x)is. m thus defined, in any of the three cases we will state 

y p (x) = [(x4 0 + A x x + . . . + A m x m )sin((3 x)+ (# 0 + + . . . + B m x m )cos(p x)|r A (25) 

k being the order of multiplicity of /(3 as a root of the characteristic polynomial associated to 
equation (21) and the coefficients {A 0 , A x , . . . , A m , B 0 , B x , . . . , B m } are the constants to be de- 
termined. 

4. If f{x) = e ax [P s {x) sin(p x) + P c (x)cos((3x)] , where P v (x) and P c (x) denote polynomials in 
x of ^ and c respectively, a, P e R . Let Mbe m = max {?,<:} , we will state 

y p (x) = e ax [(^4 0 + A x x + . . . + A m x m )sin((3 x) + (i? 0 + B x x + . . . + B m x m )cos(p x)]A (26) 

k being the order of multiplicity of a + z'P as a root of the characteristic polynomial associated 
to equation (21) and the coefficients {A 09 A l9 ...,A m ,B 0 ,B l ,...,B m } are the constants to be 
determined so that y p (x ) actually be a particular solution of equation (20). 

Let us observe that each one of these particular cases is already included in this last case. We have 
described them in this way, however, to facilitate their identification. 


Example 4.1 We want to find the complete solution of the differential equation: 
2y"{x)~ y'{x)~ y(x) = x sin(3x) + 2cos(3x) . 

For that purpose, we will first solve the homogeneous equation 

2 y”(x)-/(x)-y(x) = 0 ? 

with DERIVE, where in this case v= [2, - 1 , -1], obtaining: 


(27) 


( 28 ) 


(#1 to #33 are functions provided in CONST.MTH from DNL#33.) 


#34: DHCCCv, x) := DIFERENCIAL_HOMOGENEA_COEFICIENTES_CONSTANTES(v, x) 

x - x/2 

#35: DHCC([2, -1, -1], x) =e -c +e -c 

1,1 2,1 
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Taking into account the expression of f(x), we must find a solution of the form 

y p (*) = [ p ,n (*) sin ( 3 X ) + P m W COS ( 3 x)]x* , (29) 

being m = max{l,0} = l and k the multiplicity of3i as a root of the characteristic polynomial of equa- 
tion (28). This characteristic polynomial is given by the function AUX_1 ( [ 2 , -1 , -1]). 

2 

#36: AUX_1([2 ? -l p -1]) = 2-x - x - 1 

2 

#37: 2 ■ (3 ■ l) - 3 - 1 - 1 = -19 - 3 - i 

As we can see, 3/ is not a root of the characteristic equation obtained, so that k- 0. Thus, we look for 
a solution 

y p (x) = ( Ax + B) sin (3 x) + (Cx + Z))cos(3 x) , (30) 

where coefficients A, B, C, D have to be determined. This being a second order differential equation, 
we calculate the two first derivatives of (30): 


#38: Y_P(x) := (a-x + b)-SIN(3-x) + (ox + d)«C05(3«x) 

d 

#39: — Y_P(x) = (3-a-x + 3-b + c).C0S(3-x) - (3-cx - a + 3 - d) - SIN(3 - x) 

dx 


#40: 


f d V 

— Y_P(x) 


3 ■ (3 ■ c ■ x - 2 ■ a + 3 ■ d) ■ COS (3 ■ x) 


dxj 

Replacing in differential equation (27), we obtain 


- 3* (3 ■ a* x + 3-b + 2 ■ c) ■ SIN(3 ■ x) 


fd y d 

#41: 2- — Y_P(x) Y_P(x) - Y_P(x) 

IdxJ dx 

#42: - (x* (3 «a + 19 -c) - 12-a + 3-b + c + 19 - d) -C0S(3 -x) - (x-(19-a - 3-c) + a + 19-b + 

3 • (4- c - d))-5IN(3-x) 


Now equalling the coefficients obtained with function f(x), a system of equations is obtained, which we 
solve: 

#44: S0LVE([12 -a - 3-b - c - 19-d = 2, - 3-a - 19-c = 0, - 19-a + 3-c = 1, -a - 19-b - 

3 • (4-c - d) = 0], [a, b, c, d]) 


#45: 


a 


19 809 3 9199 

A b = A C = A d = 

370 34225 370 68450 


We check that the function obtained is a particular solution of (27) 


#46: 


Y_P(x) := 


f. _il. x _ 

809 " 

•SIN(B.x) + | 

r 3 

, v 

9199 " 

l 370 

34225 > 

# A 

, 370 

68450 > 


-C0S(3-x) 


f d y d 

#47: 2- — Y_P(x) Y_P(x) - Y_P(x) = 2-COS(3-x) + x-SIN(3-x) 

LdxJ dx 
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The general solution of (27) is given by the addition of the particular solution and the general solution 
of the homogeneous (28): 


#48: Y(x) := Y_P(x) + DHCC([2, -1, -1], x) 

x - x/2 f 3 *x 

#49: Y(x) := e «c + e *c + 

1,1 2,1 l 370 


9199 

.C0S(3-x) 

68450 J 


19- x 


370 


809 ) 

• 5IN(3 *x) 

34225 J 


Let DERIVE 6.10 do the job. 
in 0.09 seconds: 


#50: 


#51: 


DS0LYE2 



1 2 ■ COS (3 ■ x) + x- SIN(3 ■ x) 

? j 

2 2 


x - x/2 

cl«e + c2-e + 


(555 ■ x - 9199) ■ COS (3 - x) 
68450 


x 

j 

(3515. x + 161S) ■ SIN(3 ■ x) 
6S450 


[ 

© Let's try with TI-NspireCAS: 

de Solve! 2- y"-y'-y = .x m sin(3- v)+2- cosfi- xj,r\y) 

~\ 



„ 2 , x ( 3- x 9199 ' 

/ \ /-19a 809 

■ u \ 

v=c3- e +c4 ■ e + 

• cos(3- aJ+ 

■ sin (3- aJ 

\370 68450 

\ 370 34225, 



4.2 Method of variation of constants. 

The method of variation of constants or Lagrange method is more general than the method of unde- 
termined coefficients. It consists of starting from the general solution of the homogeneous problem 
associated to (20): 

y H (x)=C l y l +C 2 y 2 +... + C n y n , (31) 

and deal with the constants {Q , C 2 , . . . , C n } as functions of x, imposing them certain conditions so 
that the function thus obtained 

j // (x)=C 1 (x)yi+C 2 (x)y 2 +... + C n (x)y„ , (32) 

bo a particular solution of (20). The way to find the functions C-(x) is supplied by the following re- 
sult: 

Theorem 4.1 If [y x , y 2 , . . . , y n } constitutes a fundamental set of solutions of homogeneous differen- 
tial equation (21) associated to equation (20), and if functions {C 1 (jc),C 2 (x),... 5 C m (jc)} verify the 
system of equations: 
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C\ (x)y x + C\ {x)y 2 +... + C\ {x)y n = 0 
C'i (x)y', +C' 2 (x)y' 2 + ... + C, (x)y’ n = 0 
C\(x)y'\+C 2 {x)y" 2 + ... + C' n {x)y" n = 0 

C'| (x)y n 2) i + C 2 (x)y n 2) 2 + . . . + C\ (x)y n 2 \ = 0 

/ w 


(33) 


C, (. x)y n ~ l \ + C\ {x)y n ~ l) 2 + ... + C\ (x)y "A 

then y{x) = Cj (x)y, +C 2 (x)y 2 +... + C n (x)y n , is a particular solution of equation (20). 
Let us observe that the system of equations (33) can be written in matrix way: 


y i y 2 

y i y\ 

y" 1 y 2 


7 3 

y 3 

y 3 


y" _1) 2 y’A 


where the invertibility of the matrix 


y« 

y'n 

y"n 

„«-l) 


Y c\ (x) 2 


f 0 ^ 

c 2 (x) 


0 

c 3 (x) 

= 

0 

C (x) 


/(*) 

yv n \ A ')) 


l a n J 


(34) 


y i 

72 

7 3 

y 

y\ 

y 3 

y'\ 

y 2 

y 3 

,y n - 1} i 

y" _1) 2 

7"- 1) 


y« 

y'n 

y"n 


v »-h 
7 


(35) 


is guaranteed for {y l5 y 2 y„} being a fundamental set of solutions. In order to Apply the method 
of variation of constants, we will use the following auxiliary functions of DERIVE included in the file 

CONST.MTH 


4JY 


V_CJUJX_2(v, c, i, j) : 


ff d 

V_<LAUX_3(v, c) := VECTOR — 

UdxJ 


c, j, 0, DIM(v) - 2j 
-1 


V_0\UX_4(v, c ) := V_C_AUX_3(v, c) 

f 

V_CAUX_5(v, c, f) := VECTOR - 


■KRONECKERCj, DIM(v) - 2), j, 0, DIM(v) - 2 


V_CAU)C6(v, c, f) 
V_CAUX_7(v, c, f) 
V_CAUX_8(v, c, f) 
V_CAUX_9(v, c, f) 


V_C_AUX_4(v, c ) ■ V_C_AUX_5 ( v , c, f) 
= J V_CJVUX_6(v, c, f) dx 
= VECT0RC1, j, 0, DIM(v) - 2) 

= V_C_AUX_7(v, c, f)-c 


SOLUCION_PARTICULAR_COMPLETA( v , c, f) := V_CAUX_9(v, c, f) 
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In vector c we write a fundamental set of solutions of the homogeneous differential equation (21), 
which can be obtained by means of the function: 


DIFERENCIAL_HOMOGENEA_COEFICIENTES_CONSTANTES (v,x) . 


Function V_C_AUX_3 (v, c) calculates matrix (35). The rest of the functions serve in order to solve 
the matrix system (34). 


Example 4.2 We want to find a complete solution of the differential equation 
y m (x) + y(x) = cosec(x) . 

For that purpose, we firstly solve the homogeneous equation 
y"'{x)+ y'{x)=0 , 

with DERIVE , where in this case v = [l,0,l,0] thus obtaining the solution: 

DHCC([1, 0, 1, 0], x) = C0S(x).d + SIN(x)t + c 
#52: 2,1 2,1 1,1 

Y_H(x) := COS(x) . d + SIN(x).r + c 
#55: 2,1 2,1 1,1 

A fundamental solution set of solution of the associated homogeneous equation (37) is given by the 
terms that accompany the constants in the solution obtained: 

{l, sin(x), cos(x)}. 

We apply the function of DERIVE 

SOLUCION_PARTICULAR_COMPLETA (v, c, f ) 

where in this case v = [l,0,l,0], c = [l,sin(x),cos(x)] and /(x) = cosec(x), obtaining the particular 
solution: 

#54: SO L U CIO N_P ARTI CU L AR_CG MP L ETA ( [ 1 , 0, 1, 0], [C0S(x) ? SIN(x), 1], CSC(x)) 

#55: - LN(C0S(x) + 1) + LNCSIN(x)) - (1 - COS(x)) - x-SINtx) 

#56: Y_P(x) := - LN(CGSCx) + 1) + LN(SIN(x)) . (1 - CQS(x)) - x-SINCx) 

Adding to this particular solution the one previously obtained of the homogeneous equation, we ob- 
tain: 


(36) 

(37) 


Y(x) := - LN(CGS(x) + 1) + LN(SIN(x)). (1 - COS(x)) - x-SINtx) + CGS(x)-d + SIN(x)-r + c 
#57: 2,1 2,1 1,1 


which, as we can see, is actually a solution of equation (36): 


#53: 


Xj 


Y(x) + — Y(x) - CSC(x) = 0 
dx 
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REMARK : It is important to point out that when the dimension of the system is large, Derive will 
have plenty of problems in order to obtain the inverse of the matrix (35). This is why the method of 
undetermined coefficients is preferable. 
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In DNL#33 and #34 Rudeger Baumann posted two sequences: 

(1) [ 1 , 3 , 7 , 12 , 18 , 26 , 35 , 45 , 56 , 69 , ....] (Hofstadter-Sequence, which forms together with 

the sequence of its differences the set of the natural numbers without zero). 

and asked for a DERIVE- program to produce this sequence. 

Johann Wiesenbauer was in the mood to produce a very short DERIVE-file: 

P.S: As for the problem posed by Ruedeger Baumann in the DNL #33, here 
is a possible solution 

F (n) := REVERSE_VECTOR ( IF (n =1, [1], IF (n = 2, [3, 1], 

ITERATE (INSERT_ELEMENT (IF (SELECT (v_ = f_subl - f_sub2 + 1, v_, f_) = [], 

2f_subl - f_sub2 + 1, 2 f_subl - f_sub2 + 2), f_) , f_, [3, 1], n - 2)))) 

F ( 1 0 ) = [1, 3, 7, 12, 18, 26, 35, 45, 56, 69] 

(If you only want the nth element of this sequence, it should be easy to 
change this definition accordingly!) 


(2) The second problem was the "Look'n Say Sequence": 

[1, 11, 21, 1211, 111221, "three ones, two twos, one one" = 312211, 13112221, ] 

and in DNL#34's User Forum Rudiger wondered if it would be easy to realize this sequence in 
DERIVE, because one has to consider the order of the numbers. I tried and it took me one interesting 
evening to solve the problem (Josef). 

Look and Say: (it works, indeed!) 


H (v_) :=SELECT(u SUB l/=0 , u, APPEND (VECTOR ( IF (v_ SUB (i+l)/=v_ SUB i, [i, 

v SUB i] , [0, 0] ), i, 1, DIMENSION (v_) -1) ,[ [DIMENSION (v_) , 

v SUB DIMENSION (v_) ] ] ) ) 

H3 (v_) :=APPEND (APPEND ([ (H (v_) ) SUB 1 ], VECTOR ([ (H (v_) ) SUB (i+1) SUB 1- 
(H (v_) ) SUB i SUB 1 , (H (v_) ) SUB (i+1) SUB 2],i,l, 

DIMENSION (H (v_) ) -1) ) ) 

LOOKNSAY (st, n) : =ITERATES (H3 (u_) , u_, [st] , n) 

LOOKNSAY (1,5) = [ [ 1 ] , [ 1 , 1 ] , [ 2 , 1 ] , [ 1 , 2 , 1 , 1 ] , [ 1 , 1 , 1 , 2 , 2 , 1 ] , [ 3 , 1 , 2 , 2 , 1 , 1 ] ] 


LOOKNSAY (3, 6) 


[[3], [1,3], [1,1,1,3], [3,1,1,3], [1,3,2,1,1,3], [1,1, 1,3, 1,2, 2, 1,1, 3], [3, 1,1, 3 
,1,1,2,2,2,1,1,3] 1 


POT (v) : =VECTOR ( 10 A (DIMENSION (v) -i) , i, 1, DIMENSION (v) ) 


ZF(st,n) :=VECTOR(w . POT (w) ,w, LOOKNSAY (st, n) ) 


ZF ( 1 , 5 ) = [1,11,21,1211,111221,312211] 

ZF (1, 10) = [1,11,21,1211,111221,312211,13112221,1113213211,31131211131221, 
13211311123113112211,11131221133112132113212221] 


ZF (2, 5) 


[2,12,1112,3112,132112,1113122112] 
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The question for me - Josef - was: 

How could a function/program look like which generates the Hofstadter sequence? 

Is it possible to transfer Johann’s tricky function from above into the Nspire language? 

I found out that I would miss REVERSE, ITERATE, INSERT and SELECT among the Nspire 
functions. Reinventing these functions together with some others collected under Vector Ma- 
nipulating Functions in the DERIVE Online Help could be the first step. ITERATE will be sub- 
stituted by a loop. 

So let’s start. I create a library for manipulating lists (not vectors, because these objects can- 
not be compared with the Nspire vectors) called listops. 


Library listops containing list operations (DERIVE like): 

select1(u,k,l) returns a list of elements k of list l for which u(k) is true. 

Example: 

selectl(o>10 and isPrime(fl)=tme,fl,seq(/r,£,l J 3o)) ► { 11,13,17,19,23,29 }| 

insert(u,l,n) returns a copy of list l with expression u inserted before the n-th element. 
Example: 

insei*t(" inserted", seq (A, A, l,lo), 4 ) * { 1, 2, 3, "inserted", 4, 5, 6, 7, 8, 9, 10 } 

reverse(l) returns a copy of list l in reverse order. 

Example: 

reverse(seq(A,A,5, 15)) ► { 15,14,13,12,11,10,9,8,7,6,5 } 


select l[a> 10 and is Prime [a) = 11110 , a, seq (ft, ft, 1, 3 o)) 

{ 11, 13, 17, 19,23,29 } 

insert( 3, seq (ft, ft, 1 , 1 o), 4) { 1, 2, 3, 3, 4, 5, 6, 7, 8, 9, 1 0 } 

reverse (seq (ft, ft, 5, 1 5 ) ) 

{l5,14,13,12,ll,10,9,8,7,6,s) 

J Ei 

reverse 1/1 

Define Lib Pub re verse (/)= 
Func 

seq (/[ft], ft, dim (/), 1, - 1 ) 

End Func 

1 


selectl 6/8 

insert 1/1 

Define Lib Pub selectl (c^ >■_./_)= 
Func 

Local Hjjc 

"-(□) 

For 7 , l,dim(/_) 

/c:=c_|v_=/_[j] 

If /c=ixue 

//: =augment(//, { /_[?] }) 

End For 
// 

End Func 

Define Lib Pub insert^, ln)= 

Func 

augment(augment(seq(/[ft],ft, l,w-l),{u }),seq(/[ft],ft,w,dim(/))J 
End Func 



p 36 


The Hofstadter Sequence on TI-Nspire CX CAS 


D-N-L#35 


As there is a SELECT-function in a TI-Nspire library ( numtheory ) I took the function name 
select! Unfortunately it caused some problems which I don’t understand. The first argument 
should be a Boolean expression. The first two selectl calls are working, but combining both 
expressions by an AND fails. As you can see above one has to accomplish the second ex- 
pression to isPrime(x)=true to make it work. Don’t ask me why. numtheory\select works with- 
out this accomplishment, but only accept x as a variable and the B.e. must be entered as a 
string. 


select 1 [a > 1 0, a } seq (ft. ft, 1, 3 o) ) { 

11, 12, 13, 14,15,16, 17, 18, 19, 20, 2 1,22, 23, 24, 25,26,27,28,29,30} 

select 1 (is Prime ( z),z , seq ( ft ; ft, 1, 3 o) ) 

{2,3,5,7,11,13,17,19,23,29} 

select l(a> 10 and is Prime (x),x, seq (ft, ft, 1,3£))) 

"Error: Argument must be a Boolean expression or integer. 11 

tiitmtheo?y\se!ect(seq{k f k f l,30),"x>10 and isPrime(x) M ) 

{11,13,17,19,23,29} 

nvmtkeory\selectiseq(kX 1,30), M a>10 and isPrime(a) M ) 

"Error: Argument must be a Boolean expression or integer. " 


Now it is time to “translate” Johann’s DERIVE-threeliner. We cannot expect to do it in three 
lines, too. The IF-construct alone does need more paper work. Anf then we have to set up 
the loop (here performed as a While - EndWhile - loop). 



We compare with the DERIVE - results and can be satisfied - translation works. I wonder if 
anybody of the Nspire-programmers could do it in a more elegant way. It is not necessary to 
refer to Johann’s function but to the definition of the Hofstadter Sequence (DNL#33). Josef 



February 2001. The assembly line is buzzing and clicking. The PC department of the large 
Augsburg electronics company Compunix installs desktop and notebook computers (a) with 
Hekatium™-processors 1 2 . Supply contracts provide up to 100 processors for installation in 
each 12-hour-shift (T). The remaining material is enough to assemble 200 notebooks, a 
desktop uses 2/4 times the amount that a notebook requires (2'). A desktop computer can be 
assembled in 5 minutes, whereas a notebook computer requires 9 minutes (3'). The com- 
pany makes 150 € per desktop and 200 € per notebook (z'). 

A young businesswoman Eva Sailor is sitting in her office. The pale winter sun shines on the 
paper in front of her. She writes: 


(a) 

desktops: 

x (units) 




notebooks: 

y (units) 



(1') 

x + y < 100 

(processors) 

(1) 

y < -x + 1 00 

(2') 

y + 2.5x< 200 

(additional material) 

(2) 

y < - 2.5x + 200 

(3') 

5x + 9y < 1 2 • 60 

(assembly time in minutes) 

(3) 

y <-— x + 80 
9 




(4) 

y > 0(clear!) 




(5) 

y > 0 (ditto) 

(z') 

f t (x,y) = 1 50x + 200y 

(profits in Euro. "Maximize!" says 
the boss) 

(z) 

3 . 
y = — x + k 

4 




(maximize k) 


Ms. Sailor chews on her pencil: which production numbers are to be planned? That is which 
pair of positive integers (x,y) promises the largest profits described in the target function ( f t ) 
under the given boundary conditions (1) - (5)? 

A typical problem of the linear optimization 3 , formerly known as linear programming, a sub- 
field of operational research. This term was first used by George B. Dantzig in 1949, who 
formulated this type of problem along with its mathematical solution [3][4]. This work was 
motivated, in part, by World War II, when the U.S. Air Force required clear planning proce- 
dures for the operation of airplanes and their supplies. Naturally there had also been prelimi- 
nary studies in the civilian field, e. g. by the famous mathematician John von Neuman [1] and 
by his Russian colleague Kantorovitch [2]. In the meantime, this area of research has been 
established as a standard branch in applied mathematics. 

Our businesswoman turns round to the monitor and tries to start the company-owned PPS 4 - 
softwarepackage for production planning. Unfortunately, the in-house data network has been 
blocked for the half last hour because of heavy traffic. So, instead, she takes up her old 
school-calculator TI-92 5 and rummages through the mathematic knowledge she acquired in 
high school. 


1 This article was published first in BUS 34, Augsburg, 1998 

2 r|SKaTov (gr.) = hundred 

3 "ein Verfahren, mit dem der kleinste Oder groUte Wert einer von zwei Variablen linear abhangenden 
Funktion bestimmt wird, deren Definitionsbereich wiederum durch ein lineares Ungleichungssystem 
bestimmt ist." [6] 

4 PPS: Produktionsplanung und -steuerung [7] production planning and controlling 

5 first symbolic graphing calculator, introduced in 1995 by Texas Instruments, cf. BUS 29, july 1995, p. 

48. 
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She inputs the equations for the boundary condi- 
tions (1 ) - (4) ((5) is the y axis) in the Y= Editor 
of the TI-92 [8] (fig. 1 ). After adjusting the window 
for the coordinate system, she sees the graph 
shown in fig. 2: 

The five lines (1 ) - (5) define the region OABCD (fig. 
3). She knows that all points within or on the edge of 
this convex planning polygon are solutions (x,y) 6 
that meet the boundary conditions. She defines the 
right top border of the polygon with the function 
y5(x)=min(min(y1(x),y2(x)),y3(x)) and hatches the 
inside with the command 
Shade y4(x),y5(x),0, 150,4,6. 

Which point gives the optimum value? Surely an 
"extreme value", i. e. a point on the edge. On possi- 
bility Ms. Sailor considers to "trace" this curve y5(x) 
by typing □ Trace, and to use given coordinates to 
experimentally find an optimal point (fig. 4). 

But optimal according to which criterion? According 
to the target function (f f ) of course! She inputs the 
equation using the equation editor: y6=-.75*x+k. 

This function, however, contains the parameter k 
and therefore represents a family of lines. To look at 
the dependance of the target function of k, Ms. 

Sailor arbitrarily stores the set of values { 60 , 80 , 
10 0, 12 0} — >• k and draws the 4 selected copies 
of the target function (fig. 5). 

The desired solution is the line with the maximum 
value of k that still contains a common point with the 
set of possible solutions. As one can see, the opti- 
mal point is the vertex C. Its coordinates can be 
calculated from the intersection of lines (1 ) and (2) 
using the command 0 Math, 5: Intersection result- 
ing in the coordinates (45,55) (fig. 6). Finished. 

The procedure, however, has one weakness: The 
optimum point is determined empirically. But this 
can also be formalized. 



6 ln this case, only integer solutions are permissible, i. e. 
all grid points on or within the solution polygon 
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The polygon of possible solutions and the target 
function are mapped on the y axis as a function of 
the slope of the target function (fig. 7). Arithmeti- 
cally, a line through the origin with the slope equal 
to that of the target function is subtracted from the 
functions y4 (lower edge) and y5 (top margin of the 
planning polygon) (in the # Editor the equations y7 
and y8). By this geometric mapping the target func- 
tion becomes horizontal (that is the trick!). The op- 
timum point becomes the maximum of the path y7. 
Its maximum C can easily be determined automati- 
cally by F5 Math, 4: Maximum (fig. 8). 

Mapping back C on the original planning polygon to 
C is only a formality. Its coordinates (x,y), fortu- 
nately integers, now supply the optimal production 
numbers for one shift: 45 desktop computers and 55 
notebook computers and a maximized profit of 
z(45,55) = 17750€(fig. 9). 






^ 5_ 



n 

1 

Iff 

1A .. 

xc: 4, 

>• 

\ yc:55. 

LINDF'TIM 

HD flUTD FUNC 


Our businesswoman puts her TI-92 aside and leans back. The whole procedure took hardly 
more than 8 minutes. She solved the linear optimization problem graphically - using a nu- 
meric tool. The problem is representative for a set of similar problems. Their solutions can be 
integers as in this case or real numbers. The problems are usually subject to more than five 
boundary conditions and lead in general to solution polygons with correspondingly more ver- 
tices. The target function may need to be maximized or minimized and the solution may be a 
single point (vertex) or a whole segment (edge). 

Further examples include transportation problems (the famous "problem of the travelling 
salesman"), menu planning, traffic planning, controlling of machine production, and above all, 
economic optimization. In all cases the occurring relations are either linear or can be ap- 
proximated as linear. 


The graphic solution presented here can only be applied in this form with two variables (x 
and y). In the case of three variables, the diagram is already more complex: 


consider the target 
function: 


(z) f z (x,y,z) = x + 1 0y + 1 0Oz 


with the boundary 
conditions: 


(1) 4x + 25y + 5z< 100 

(2) 25x + 5y + 4z < 1 00 

(3) 1 0x + 3y - 30z > - 30 

(4) 25x + 8y + 1 0Oz < 200 



(fig. 10) 
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These inequatilities describe semi-spaces, which 
are cut off by planes in three-dimensional space. 

The solution space, given by the intersection of the 
semi-spaces, is a polyhedron (fig. 11). 

Even here, the TI-92 can be used for visualization, 
although it scarcely offers a numeric assistance 
during the optimization with the graphic procedure. 

Using F3: Trace, one can trace on the grid lines of 
the polyhedron and approach the solution, deter- 
mined by the resolution of the grid (fig. 12). How- 
ever, the identification of the correct corner point is 
left to the eye. 

A trick helps here also. If one rotates the axes such 
that the plane of the target function is horizontal, 
then the maximum shows up as "the highest eleva- 
tion of the solution mountains" (fig. 13). 

Inequalities with more than 3 variables cannot yet 
be visualized, not even in the year 2001 . They rep- 
resent hyperpolyhedrons in a space of the dimen- 
sion higher than 3, which are limited by hyperlevels 
of one dimension less. For such optimization tasks 
one remains set back on the simplex procedure [5], 
a numeric algorithm which goes back to George B. 

Dantzig [4]. 

Then the telephone rings. It’s the department manager. "Ms. Sailor, I just remembered we 
have two assembly lines that are no longer being used. What if we convert them to ...? Can 
you calculate it again? ... - say max. 15 per shift - and above that the processors for the new 
pocket calculators... ". With a soft "ping" Ms. Sailor’s monitor reports back. In the menu of her 
PPS software she selects... linear optimization... simplex procedures... 
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I append the graphic approach for treating 2 dimensional linear programming problems with 
the Tl-Nspire CX CAS: 
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The procedure is very similar to the way to work with the TI-92. Unfortunately I was not able 
to reproduce the 3D-representation of Detlev’s second problem. Josef 


□ f^)=\ 
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Die Platonischen Korper - Platonian Solids - 
- Hidden Lines with TI-92 - 

by Heinz Rainer Geyer, St. Katharinen, Germany 
HeinzRainer . Gey er @t-online . de 

Inspired by Hubert Wellers "s presentation of poly- 
hedrons (see DNL #29, p20-p25) we spent a lot of time 
in our school (in the scope of a working group CoMa - 
Computer-Mathematics - at the Gutenbergschule in 
Wiesbaden) with Platonian Solids. 

A dodecahedron in isometric representation > 

It was easy to realize the Hidden-Line- Algorithm using DERIVE. Working with the TI-92 we had to 
overcome some ’’hurdles”, because the needed matrix operations - better operations applied on lists of 
matrices - are not possible. Creation of a new matrix from rows of an existing one has to be done ele- 
ment by element. Fortunately all boundary planes of Platonian Solids show the same number of verti- 
ces, so we could easily solve this problem. 

We begin defining the solid by one matrix of 3D-coordinates. This matrix contains the closed se- 
quence of points of each boundary plane. Matrix icosa(a) for the icosahedron consists of 20*4 = 
= 80 rows (20 planes with (3+1) points each). 

In our realisation variable a has various meanings: it is the length of an edge (Tetrahedron and cube) 
or the length of the base-cube's edge (octahedron, dodecahedron and icosahedron). 

Entering the matrices based on our DERIVE- results was a bit laborious but it was also a useful training 
for our concentration and finally we obtained the matrix. I will not explain now the calculation of the 
coordinates. 

The program allin(body,proj,n) draws all edges 
of solid body using the matrix of projection proj. 
n is the number of rows per face. 

(Try: allin(dodeka(4),kav1 ,6)) 

Matrix proj converts matrix body by multiplication to 
a matrix representing a 2-dimensional object - a set of 
points in the plane -, which are all drawn using the line - 

- command. 

hidd(body,proj,n) performs for each face the determi- 
nant-test before drawing. Because of the problems men- 
tioned above the determinants are calculated in a subrou- 
tine. ( hidd(dodeka(4), kavl , 6) ) 

(Don't forget: the sequence of points defining the faces 
need positive orientation). 



FEt 

fFsy 


r fe 

T FfiT T 'i 

Vi — p 

untrol 

I/O 

Oar 

Find 

. . |Mode| | 

: hidckr 

pro 

n J 




■ Prgn 
: Local 

n, 1 , i . 

, tes 

■t 



:0 Tesl 

.routine u 

iith 

_Dete 

rn inant 


: Define testtmm, k>=Func 
: tnn [ k+2, 1 ]-mm [ k+1 - 1 ] [ k 7 2 ]-rm [ V 
])-<m[t+2 J 2]-m[t+l,2]>*<nn[k, 1 1-p 
1 , 1 ] > 

: EndFunc 

:0 Start MainProg 

!(i Trancfoi" TH— P.-. i h-.+ ■= + VT I 

: n3*pro->n 

:For 1 ? 1 ? rowDin<rO ? n 
: If testt 1 ) >0 Then 
: For i ? G, n-2, 1 

: Line pitl + i, 1 ],pi[l + i,2],m[l + i + l, 1 ] 

+ i + l 7 2 ] ? 1 
: EndFor 

: End If 
: EndFor 
: EndPrgrJ 




5 DLID 5 


RAD AUTD 


FUNC 
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On the defined solids you can apply appropriate projections. Multiplications of matrices can be per- 
formed without any problems. 

If you want to translate the solid you should add a constant vector v of same dimension to each row 
body[i] of the matrix. On the TI-92 PLUS such program works, the ordinary TI-92 produces a 
Domain Error, although a singular calculation body[i]+v can be performed on the Home-Screen. 
Again the solution is an addition element per element which can easily be followed in the program 
trans(body,vec): The solid body is translated by the vector v. 


rFiTmf fit yfstVfhty fs 

|t IControl |l/0|UarjFind. , 

T T i 

. . |Mode| 1 

: axestfcn, yn ? zn ? pro) 


: Prgn 


: Local n 


: [ [ xn , 0 , 0 ] [0,0,01 [ 0 , yn , 0 ] [0,0,01 [ 0 , 0 , zn ] 

]*pro*n 


: plotn(n) 


: EndPrgn 



IS 


YFStYfHtY FE 

Control I/O Oar Find. 


trans < n , o ;■ 

Func 
Local i ? k 

For i ? 1 ? rowDin^n) ? 1 
For k, 1 ? colDin<o) , 1 
n [ i , k ]+o [ 1 ? k ]->n [ i ? k ] 
EndFor 
EndFor 
Return n 
EndFunc 


SOLIDS FiliD ftUTO FUNC 


Program axes(xm,ym 5 zm 5 pro) represents the axes. Parameters are the oriented lengths of the three 
axes and the projection matrix in use. 




Ikosaeder in ’’Military Perspective” Cube together with a translated 

Octahedron in dimetric projection 




|Rlgebra|Calc|Other|FrgpiIO|Clear a- 


■ axest? ? 7 ? 7 , ni 12) 

Done 

■ hidd<; icosa<;5) , ni 12 , 4) 

□one 

■ axes<;7 ? 7 ? 7 ? dine) 

□one 

■ hidd<;cube(5) 7 dine 7 5) 

□one 

■ hiddi;trans(;okta(5) ? [2. 5 

0 0]) ? dine ? 4) 


Done 


Using all these tools one can do a lot in Linear Algebra. 
Working with nonregular solids you will miss DERIVE r s 
wonderful VECTOR-command. 


The projection matrices used here are: 



— 1/2 

-1/2' 


'S/2 

-1/2 ’ 


~S/2 

-1/2 " 

kavl = 

1 

0 

; mill = 

1/2 

S/2 

; mil2 = 

1/2 

S/2 


0 

1 


0 

1/5 


0 

1/2 
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Heinz Rainer wrote about nonregular solids. I faced that challenge and modified - extended - his pro- 
grams for presenting other solids on the Tl-screen. 

Let's take a pyramide with a hexagonal base. 

The base is defined by six points (A,B,C,D,E,F). We first enter each single point as a list consisting of 
its three coordinates, followed by the vertex S(0,0,4): 


f'Fi T^Y FET Y FjT Y FHt Y FE Y Ffi T T "l 

E |fllgebra|Calc|0ther|PrgmI0|ciean Up| 1 

■ £4 

G -2J + a : 

£2 

2 13 -2Hb i 

£► 




£ -2 

2-J3 

-2> 

■ £ -4 

G -2>+d 

: £ 

2 "2 f3 

-2> - 

e ► 




£2 

-2-.J3 

-2> 

■ £0 

G 5> + s 



£Q G 

5> 

create<31> 1 

GEYER 

RAD EKACT 

FUNC 3/30 



ifiT^T FT? fht y— 

f— |Rlgebra|Calc|Other|FrgroIO|Clean 


ifc lean Up^l 


Name of solid: Ihexpyrl" 





, £Enter=0K 3 

£ESC=CFlNCEL!) ^ 



£2 -2-TS -2> 

■ £G G 5> s 

£Q G 5> 

create<31> 1 

GEYER RAD EKACT 

FUNC 3/30 


As Heinz Rainer mentioned before editing a matrix consisting of 7 + 6*4 = 31 points is very boring. I 
use the program create ( number_of_poi nt s ) to produce the matrix, describing the pyramid: 


Take care of the order of the points and enter 3 1 points: 


A,F,E,D,C,B,A, A,B,S,A, B,C,S,B, C,D,S,C, D,E,S,D, E,F,S,E, F,A,S,F 





. FE YF3 , 

iPloi SetupjCel 1 1 


MAT 

31x3 





■ 

cl 

c2 

c3 

c4 


1 

4 

Q 

-2 



2 

2 


-2 




-2 


-2 




-4 

G 

-2 




-2 

2*T<3) 

-2 




2 

2*r<3) 

-2 





Q 

-2 


■ 


rlcl=4 


If you want to change some values, then open matrix 
hexpyr with the Matrixeditor and perform your 
changes. Define a vertex- vector with the numbers of the 
faces' vertices and then call 

al 1 i ngen (body_name , proj ecti on , 
v_vector) or 

hi ddgen (body_name , proj ecti on , 
v_vector) . 


fFlTTlT FEt Y F3t Y FHt Y FE Y Fb^ 

| t f— |fl 1 gebra |Ca 1 c |0ther |Prgm 1 0 |C 1 ean 


■£^4 G - 2> + d : £-2 +TYT 

-2> +e : 

£2 

-2 J3 

-2> 

■ £G G 5> * s 

£Q 

G 5> 

■ create£31) 


Done 

■[6 3 3 3 3 3 3] + faces 

[6 3 3 

3 3 

3 3] 

■all ingeKhexpyr ? kaol ? faces) 


□one 

■ h iddgenC hexpyr ? kaul ? faces) 


Done 

hiddsren Oiexpyr , kaul , f aces > 

GEYER RAD EKACT FUNC ?/30 




Try some other projections: 



'-1 -1/2" 


i 

0 

iso = 

1 -1/2 

; or oblique view obl(sf, O) = 

sf cos( O) 

-sf sin( O) 


0 1 


0 

1 
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The Tl-programs: 

create(num) 

Prgm 

Local j_, num , bo , pt 

Request "Name of solid",nams 

{}-bo 

ClrlO 

For j , 1 , num 

InputStr "Point " , pt 
augment (bo, expr(pt) )^bo 

EndFor 

1 i st >-mat (bo , 3)^bo 
bo^#nams 

Disp "ready, press F5" 
EndPrgm 


hiddgen (m3 , pro , f ) 

Prgm 
Cl rDraw 

Local m , 1 , i , j , st , test , k_, f 1 a 
© Testroutine with Determinant 
Define test(mm,k)=Func 

(mm[k+2 , 1 ] -mm[k+1 , 1 ] ) * (mm[k, 2] -mm[k+1 , 2] ) - (mm [k+2 , 2] -mm [k+1 , 2] ) * (mm [k, 1 ] - 
mm[k+1 , 1 ] ) 

EndFunc 

ClrlO 

© Start MainProg 
© Transfer 3D-Points to 2D 
m3*pro^m : l^st : dim(f ) [2]->f 1 a 

For k , 1 , f 1 a , 1 

If test(m,st)>0 Then 
For i , 0 , f [ 1 , k_] -1,1 

Line m[st+i , 1 ] ,m[st+i ,2] ,m[st+i+1 , 1 ] ,m[st+i+1 ,2] , 1 
EndFor 
Endlf 

st+f [1 , k_] +1^st 

EndFor 

EndPrgm 

For allingen(m3,pro,f) you have only to remove the hidden line test. 

I add another object stairs (produced with a vector fstair). 
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On many occasions I am asked how to produce 3D-pictures of solids of revolution on the TTs 3D- 
screen. Students often show an enormous deficiency in their 3D-imagination. Producing their own 
representations might improve their visualization of 3D-solids. The procedure is very easy to perform. 
A simple sketch describes how to do it. 

Let's take any function v - fit) which should be rotated around the v-axis. Then the v-axis changes to 
the z-axes and argument t = ^ fx 2 + y 2 . Function z(x, y) = f (^J x 2 + y 2 ) describes the solid. 


For an easy example we want to rotate the parabola y = — around the y-axis. 


✓ 



1. Step: 

In the Home Screen define function 

x A 2/2 k(x). 

2. Step: Switch the Graph Mode to 3D. 

3. Step: 

Define in the #-Editor the auxiliary function zl (x,y) 
and then the surface z2(x,y). Set the 3-values. Press 
% and wait for the graph. 





The next example is more difficult: we take a 
piecewise defined function and want to produce 
the solids resulting from rotations around both 
axes. 


The function is composed of two parabolas: 
Vx for 0 < x < 4 

y = < 4- ^/8 - x for 4 < x < 8 
undefined else 



when(x<0, undef,when (x < 4, V (x) , when (x <8,4 - V(8-x) , undef ) ) ) 





p 48 


Josef Bohm: Solids of Revolution on the TI-92/89 


D-N-L#35 






iRLDTS 

z2=k(zl(x ? y)) 
^ z3= /u ndef ? zU x, y)< G 


t -Jzl(x ? y),zl(x, y ) <4 . 

f4 - JS - zl(x 7 y) ? zl(x 7 y) ^ 3 , ? 

\_mdef ? else ?ei£e 

1 undef,zl(x, 

f(zlCx ? y)) 2 ; 


■ y)< 0 
,zl(x, y)< 2 


z3<x J .^>=nhen<al<x J .k><D J .iindef 

MHIH FiliD HUTD 3P 




Rotating the curve around the x-axis should result in a bell- 
shaped surface. 

In the #-Editor edit function z3(x,y) according to the func- 
tion given above , replacing x by zl(x,y). 

Take care of the 3-settings and plot the graph. 


As this is not a function, we must use a trick: we rotate the inverse function around the y-axis to obtain a 
visualization of the bell. We also have to adapt the 3-settings. Take into account that the scaling is not the 
same in both windows. So the bell doesn't look as slim as it should. 






"4 - JS - zl(x 7 y) ? zl(x 7 y) ^ S , 
.mdef,else ?ei£e 


, z4 _|undef ? zl(x ? y)< G 
" 1 fCzl(x ? y)) 2 ? zl(x ? y)< 2 


z5= 


"(zlCx, y!VP + 0-zlCx, y) - 0,zl(x , y 
_undef , else 


z4<x J .y)=,,,8*zl<x J .y>-8 J .imdef>>> 

MHIN RHP HUTD 5P 





MHIN RHP HUTD 5P 


Finally let me show, how the "Chinese Hat" and the "Bell" appear in DERIVE. 



The right picture is produced with an a- version of DERIVE 5, which was presented by David Stoutemyer 
and Theresa Shelby at the ACDCA-Summer Academy. It should be an appetizer for the new DERIVE- 
age. 
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Let’s try to present the solids of revo- 
lution on the TI-Nspire CX CAS 
screen. The rotation of the parabola is 
an easy job - the procedure performed 
on the TI-92/Voyage 200 can be trans- 
ferred on the Nspire without any prob- 
lems: 




I was not able to repeat the procedure 
with the piecewise defined function - 
the Chinese Hat and the Bell. 


I failed with all my tries and efforts. 


If anybody of you will have more suc- 
cess, then please keep me informed. 


*w-- 

Fertig 

4 



Fertig 

z2{?,y)=k(zl(x,y)) 

Fertig 


xV 


4 

J5T. o<x<4 

Fertig 

kl(x}=U-ft^ 4 <x<8 


(undef, x<0orx>8 


z3(xj/):«*7(z;M) 

Fertig 

p7(x):=^x 

Fertig 

T3 

£ 

ll 

tl 

Fertig 

PM-* 2 

Fertig 

pM- 8-(4-x) 2 

Fertig 


Next experiment was to define the two parts of both solids in parameter form. This works as you can see 
on the next screen shot. 
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pi (t) and p2 (t) form the red solid 
and p3 (t) with p4 (t) form the blue 
one. 



Using an idea of Michel Beaudin 
and his colleagues from Montreal 
we can “invent” DERIVE ' s CHI- 
function for defining a piecewise 
defined function. 

Then we can use this function 
which is called kl(x) in my screen 
to plot the z2(x,y ) created solid 
(rotation around the z-axis) and the 
parameter form created solid (rota- 
tion around the x-axis) using only 
one function. 




